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DISCRETIZATION OF DIV-CURL SYSTEMS BY WEAK GALERKIN 
FINITE ELEMENT METHODS ON POLYHEDRAL PARTITIONS 

CHUNMEI WANG* AND JUMPING WANGt 


Abstract. In this paper, the authors devise a new discretization scheme for div-curl systems 
defined in connected domains with heterogeneous media by using the weak Galerkin finite element 
method. Two types of boundary value problems are considered in the algorithm development: (1) 
normal boundary condition, and (2) tangential boundary condition. A new variational formulation 
is developed for the normal boundary value problem by using the Helmholtz decomposition which 
avoids the computation of functions in the harmonic fields. Both boundary value problems are 
reduced to a general saddle-point problem involving the curl and divergence operators, for which 
the weak Galerkin finite element method is devised and analyzed. The novelty of the technique lies 
in the discretization of the divergence operator applied to vector fields with heterogeneous media. 
Error estimates of optimal order are established for the corresponding finite element approximations 
in various discrete Sobolev norms. 

Key words, weak Galerkin, finite element methods, Helmholtz decomposition, weak divergence, 
weak curl, div-curl systems. 
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1. Introduction. This paper is concerned with new developments of numerical 
methods for div-curl systems with two types of boundary conditions. The model 
problem seeks an unknown function u = u(x) satisfying 

(1.1) V • (piu) = /, in n, 

(1.2) V X u = g, in 12, 

where 12 is an open bounded and connected domain in with a Lipschitz continuous 
boundary T = 912. The Lebesgue-integrable real-valued function / = /(x) and a 
vector field g = g(x) are given in the domain 12. Here fi = {/iij(x)} 3 x 3 is a symmetric 
matrix, uniformly positive definite in 12 and with entries in L°“(12). Assume that the 
boundary T has m -I- 1 connected components such that 

m 

r = U r.. 

i=0 

where Tp represents the exterior boundary of 12, and Ti,l < i < m, the other con¬ 
nected components of F. 

The system (fFTt-lfL^ arises in fluid mechanics and electromagnetic field theories. 
In the fluid mechanics field theory, the coefficient matrix ^(x) is diagonal where 
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diagonal entries are the local mass density. In electrostatics field theory, /i(x) is 
the permittivity matrix. In linear magnetic field theory, the function /(x) is zero, 
u represents the magnetic field intensity and /r(x) is the inverse of the magnetic 
permeability tensor. 

We consider two types of boundary conditions for the div-curl system (fra-lfOl): 
the normal boundary condition, and the tangential boundary condition. 

The normal boundary condition is concerned with a given flux value for the vector 
field on T; i.e., 

(1.3) inu) • n = ^, on r, 

where n is the unit outward normal direction on T. 

The tangential boundary condition corresponds to a given value for the tangential 
component of the vector field u; i.e.. 


(1.4) 


u X n = X, on r, 
(^u • rii, l)ri =/3i, i = 


where is the unit outward normal direction on the connected component T^. 
The space of harmonic fields is given by 

H^(II) = {v e : Vxv = 0, V • (^iv) =0, /rv n = 0 on T}. 


The space of harmonic fields is non-trivial when the domain is not simply connected. It 
is readily seen that the div-curl system (HID-diiD with the normal boundary condition 
(lOl) is generally not well-possed, as uniqueness fails (adding a harmonic field to a 
solution u still makes a solution). For this reason, we shall look for a solution u that 
is orthogonal to the space of harmonic fields in the /r-weighted norm. Specifically, 
the complete problem reads: given g € [L^(II)]^ with V ■ g = 0 in fl, find u such that 


(1.5) 


' V ■ ifiu) = /, 
V X u = g. 




/ru • n = 5, 
/ru • J 7 = 0, 


in fl, 
in n, 
on r, 

yrj 


The problem (II.5p has a solution [55], and uniqueness is straightforward. It is also 
well-known that the div-curl system (frT])-(lT2l) with the tangential boundary condition 
dLl has one and only one solution. 

There have been many numerical methods for approximating div-curl systems. In 
[22] . Nicolaides proposed and analyzed a control volume method directly for planar 
div-curl problems. In [53], Nicolaides and Wu presented a special co-finite volume 
method for div-curl problems in three dimension, which was based on a system of two 
orthogonal grids like the classical Voronoi-Delaunay mesh pair. Bossavit [3] proposed 
a classical numerical method for solving the magnetostatic problem by introducing 
a scalar or vector potential. In |12j . a discrete duality finite volume method was 
presented for div-curl problems on almost arbitrary polygonal meshes. In [4], Bramble 
and Pasciak proposed and analyzed a direct numerical scheme under a very weak 
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formulation where the solution space was In [TU], a mixed finite element 

method was analyzed for div-curl systems in simply connected axisymmetric domains. 
In [26] , a numerical algorithm was designed for constructing a finite element basis for 
the first de Rham cohomology group of the computational domain, which was further 
used for a numerical approximation of the magnetostatic problem. In [7] and 
the mimetic finite difference method was applied to the 3D magnetostatic problems 
on general polyhedral meshes. 

Recently, weak Galerkin (WG) finite element methods have emerged as a new 
numerical technique for approximating the solutions of partial differential equations. 
The WG method was first introduced in [^[32] for the second order elliptic problem 
and was further developed in [sni na EH [ 29 ] with other applications. Two basic 
principles for the WG finite element method are: (1) the differential operators (e.g., 
gradient, Laplacian, Hessian, curl, divergence etc.) are interpreted and approximated 
as distributions over a set of generalized functions, and (2) proper stabilizations are 
employed to enforce necessary weak continuities for approximating functions in the 
correct topology. It has been demonstrated that the WG method is highly flexible 
and robust as a numerical tool that makes use of discontinuous piecewise polynomials 
on polygonal or polyhedral finite element partitions. 

The goal of this paper is to present a new discretization scheme for the div- 
curl system (iTT])-(fT:^ in any connected domain with heterogeneous media by using 
the weak Galerkin finite element approach. In particular, for the normal boundary 
value problem (USD, a new variational formulation is developed by using a Helmholtz 
decomposition which avoids the computation of any harmonic fields v S 11^(17). 
The resulting formulation is a special case of the model problem (15.11) detailed in 
Section 5. The div-curl system with the tangential boundary condition dm is also 
formulated as a special case of the model problem dm- Therefore, our attention is 
focused on the development of weak Galerkin finite element methods for (15.111 . It is 
readily seen that the main difficulty in numerical methods for dSD) lies in the term 
V • (/ru) which requires the continuity of fj,u in the normal direction of any interface, 
particularly the interface of any two polyhedral elements. The weak Galerkin finite 
element method offers an ideal solution, as the continuity can be relaxed by a weak 
continuity implemented through a carefully chosen stabilizer. 

The paper is organized as follows. In Section 2, we introduce some commonly used 
notations. In Section 3, we derive a formulation for the div-curl problem dni)-dii2]) 
with normal boundary condition dm by using Helmholtz decomposition. Section 4 
is devoted to a discussion of the div-curl system with tangential boundary condition. 
In Section 5, we discuss a model problem that is central to the solution of the div- 
curl system with both the normal and tangential boundary conditions. In Section 
6, we introduce some discrete weak differential operators which are necessary for 
the development of weak Galerkin finite element methods in Section 7. Section 8 
is devoted to a discussion of existence and uniqueness for the solution of the weak 
Galerkin discretizations. In Section 9, we derive some error equations. An inf-sup 
condition is established in Section 10. Finally in Section 11, we derive some optimal 
order error estimates for the WG finite element approximations. 

2. Notations and Preliminaries. Throughout the paper, we will follow the 
usual notation for Sobolev spaces and norms [5]. For any open bounded domain 
D C with Lipschitz continuous boundary, we use || • ||s^£) and | • to denote the 
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norm and seminorm in the Sobolev space H‘^{D) for any s > 0, respectively. The inner 
product in H’^[D) is denoted by (•, ■)s,d- The space H^{D) coincides with L‘^{D), for 
which the norm and the inner product are denoted by || • lie and (•: ■)d, respectively. 

Let ^ = {^iij} 3 x 3 be a symmetric matrix, uniformly positive definite in D and 
with entries in L°°{D). Introduce the following Sobolev space 

iL(div^; Z?) = {v G [L\D)f : V ■ (/iv) G L\D)}, 

with norm given by 

l|v||«(div,;e) = (l|v||!, + ||V.(/iv)|||,)^ 

where V ■ (/rv) is the divergence of /rv. Any v G iZ(div^; D) can be assigned a trace 
for the normal component of /rv on the boundary. The subspace with vanishing trace 
in the normal component is denoted by 

iZo(div^;Z:>) = {v G iZ(div^;II) : (^v) • n|a£, = 0}. 

Denote the subspace of iZo(div^; D) with divergence-free vectors by 

F^(D) = {v G iZ(div^; D) ; V ■ (/rv) = 0}. 

When fj. = I is the identity matrix, the spaces iZ(div^;D), iZo(div^;D), and F^(D) 
shall be denoted as iZ(div;D), iZo(div;D), and F(D), respectively. 

Denote by iZ(curl; D) the following Sobolev space 

iZ(curl;D) = {v : v G [L^{D)f,V x v G [L^{D)]^} 

with norm defined by 

l|v||^(curl;D)=(||v|||, + ||Vxv|||,)A 

where V x v is the curl of v. Any v G ZZ(curl; D) can be assigned a trace for its 
tangential component on the boundary. The subspace of ZZ(curl;D) with vanishing 
trace in the tangential component is denoted by 

ZZo(curl;D) = {v G ZZ(curl;ZZ) : v X n\gD = 0}. 

When ZZ = D, we shall drop the subscript D in the norm and inner product 
notation. 

For simplicity of notation, throughout the paper, we use “< ” to denote “less 
than or equal to up to a general constant independent of the mesh size or functions 
appearing in the inequality”. 

3. The div-curl System with Normal Boundary Condition. The goal 
of this section is to derive a suitable variational formulation for the problem dun. 
Denote by ZZ°(curl;D) = {v G : Vxv = 0} the space of curl-free fields. It 

is well-known that any vector field v G ZZ°(curl; ft) can be written as (see, e.g., [27) 1 

(3.1) v = \7(j) + ri, 

where </> G ZZ^(D) and r) G IH^(D). It follows that 

(V((), fiT]) = -{(j), V ■ (/rr;)) -h (0, /rr; • n)r = 0. 

The decomposition (EH) is thus orthogonal in the /r-weighted norm in ZZ*^(curl; D). 
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3.1. Helmholtz decomposition. Denote by the set of functions in 

with vanishing value on Fq and constant values on other connected components 
of the boundary; i.e., 

= {'/'€ H^{n) : (j)\ro = 0, = c, i = 

Next, introduce a Sobolev space 

Y^(D) = {v G iFo(curl; D) : (/rv • n^, l)ri =0, i = 1,. .., m}. 

The following Helmholtz decomposition holds the key to the derivation of a suitable 
variational form for the problem (11.51) . 

Theorem 3.1. For any vector-valued function u G there exists a unique 

if G Y^(D) nF^(r2), (f G iF^(r2)/R, and r] G H^(D) such that 

(3.2) tt =X i/j + + r/. 

Moreover, the following estimate holds true 

(3.3) ||’MlH(curl;n) + l|V(()||o < {kU, u)^ . 

Proof Consider the problem of seeking xf G Y^(H) D F^(H) such that 

(3.4) (/r-^V X 7/),V X <^) = (u,V X ¥>), V ¥> G Y;,(H) n F;,(H). 

Denote by 

a{'ip,ip) := (^“^V X xp,S/ X ip) 

the bilinear form defined on Y^(H) D F^(H). We claim that a(-, •) is coercive with 
respect to the iF(curl; D)-norm. To this end, it suffices to derive the following estimate 

(3.5) l|v||o< llVxvIlo, VvG Y^(D)nF^(H). 

In fact, for any v G Y^(D) D F^(D), from Theorem 3.4 of |13) . there exists a vector 
potential function uj G [7F^(D)]^ such that 

(3.6) /rv = V X Ci), V • a; = 0, ||u;||i < (/rv, v)^. 

Using the integration by parts and the condition v x n = 0 on F, we have 

(/rv, v) = (V X u:, v) = (a;, V x v). 

It follows from the Cauchy-Schwarz inequality and (13.61) that 

(l‘v,v) < Iloilo ||V X v||o < (/iv,v)5 ||V X v||o, 
which implies (13.51) . 

Now from the Lax-Milgram Theorem, there exists a unique xp G Y^(f2) D F^(r2) 
satisfying the equation p.4l) such that 

||V’l|H(curl;n) < ||u||o < (/rU,u). 
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It is easy to see that Y^(0) n F^(0) is equivalent to the following quotient space: 

i?o(curl; n)/(Vi?o^(ri)) = {v e i7o(curl; fl) : (^v, V(j)) =0, V(/) G Hq^{Q)}. 

Thus, by using a Lagrange multiplier p G HQ^{n), the problem (13.41) can be re¬ 
formulated as follows: Find xp G iLo(curl; fl) and p G iLoc(^) such that 

X ■0, V X <^) + (/rVp, ¥>) = (u, V X v?), V v? G iLo(curl; fl), 
( 0 ,AiVs)=O, VsGiLi,(r!). 

It follows from the first equation of (13.7p that 

V X (u — X 0) — pS/p = 0. 

Since p G then the two terms on the left-hand side of the above equation are 

orthogonal in the /i“^-weighted norm. Thus, 

V X (u - X 0) = 0 ^ u - X 0 G iL°(curl; fl). 

Thus, there exist unique (p G 7L^(r2)/M and r/ G H^(ri) such that 

u — K~^y X 0 = yp -|- 77 , 

which completes the proof of the theorem. □ 

3.2. A variational formulation. Assume that the div-curl problem (inii-dia 
with boundary condition m has a solution. Integrating (HID over the domain ft 
and from the integration by parts, we have 

if, 1) = (V • (pu), 1) = {{pu) ■ n, I)r. 

Using the boundary condition m, we arrive at the following compatibility condition 

(3.8) (/,l) = (e,l)r. 

In addition, taking the divergence to the equation (HID, we obtain a second compat¬ 
ibility condition 

(3.9) V.g = 0 inO. 

Lemma 3.2. Let u be a solution of the div-curl system JO)-JO) with the bound¬ 
ary condition GJ). Then, u can be decomposed as follows 

(3.10) 7 t =X 0 -I--I- 77 , 

where 0 G Y^(r2) (~l F^(r2) is the unique solution of the following equation 

(3.11) (M-'Vx0,VxC) = (g,C), VCe Y^(f^)nF^(U), 

and (j) G satisfies 

(3.12) {pyp,yv) = {^,v) - {f,v), yv€H\n)/R. 
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Proof. Using the Helmholtz decomposition dH in Theorem o there exist 
unique 'll] G Y^(r2) fl F^(r2), cj) G and r] G H^(H) such that 

(3.13) u = X i/) + V()) + ry. 

For any C S Y^(H) n 

(V(/), VxC) = (VxV<)), C) - (V0, C X n)r = 0 . 

Thus, by testing both sides of (13.131) with Vx^, we obtain 

(3.14) (u - rj, VxC) = ip-^\7xxp,yxC). 

From the equation (HU and the fact that Vxr] = 0 and C x n = 0 on F, we obtain 

(u-r 7 ,VxC) = (Vxu,C) = (g,C)- 
Substituting the above into (13.141) yields 

(3.15) (M-'Vx^,VxC) = (g,C), VCGY^(H)nF^(H), 
which verihes ()3.11l) . 

Next, we test (13.131) against any fiS/ip to obtain 

(3.16) (u, ^V(/5) = (^“^Vx'j/),/xV(/3) + (/rV(/), V(/j) + (j 7 ,/rVvj), G i/^(H). 

Using the integration by parts, 

(u, fiVip) = -(V • (mu), (^) + ((mu) • n, (p)r. 

Thus, from the equation CH) and the boundary condition m, we have 

(3.17) (u,mV(m) =-(f,<M) + (^,(M)r. 

Similarly, from the integration by parts and the fact that rj G H^(H), 

(3.18) (i 7 ,mV^) = -(V • (Mr7),<M) + {(m^) •n,(M)r = 0. 

Since if & Hq{cut\; H), then 

(3.19) (Vx-i/), V(m) = 0. 

Substituting p. 17^ - 03.19p into p.lGp gives rise to 

(3.20) V(m) = (e, ‘M)r - (f, ‘m), V(m G H^^). 

The above problem has a unique solution in the quotient space i7^(H)/R due to the 
compatibility condition (13.81) . This completes the proof. □ 

By using a Lagrange multiplier p, the problem ()3.11|) can be re-formulated as a 
saddle point problem that seeks if G Y^(H) D H) and p G satisfying 

(m-'VxV^, VxC) + (V • (mC),p) = (g, C), 

(V ■ {pif),w) = 0, 
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(3.21) 


VCe Y^(H)ni7(div^;H), 
V u; G L^{n). 


















Going back to the well-posed problem dm, it is easily seen that the solution of 
(11.51) also admits the decomposition (13.101) . Thus, using ip x n = 0 and fj,T] ■ n = 0 on 


r, we obtain 


0 = 

by last condition in (11.51) 

= (VxV^, rj) -b (mVC), rj) -f (m??, rj) 

bv the decomposition (13-lOD 

= Vxrj) - (Cl, V • (mi?)) + (m’7,'^) 

by integration by parts 

= (mi?,^)- 

by ?7 e H^(fl) 


It follows that t] = 0. The result can be summarized as follows. 

Theorem 3.3. Let u be the solution of the div-curl system Then, u can 

be represented 

(3.22) tt = V x^p + V(p, 

where xp G Y^{Ll) n F^(r2) is the unique solution of the system of equations \3.21\) . 
and (p G is determined by the equation \3.1‘A) . 

Theorem l3.3l indicates that a suitable variational formation for the div-curl system 
dESl) is given by (I3.2ip and (13.121) , which are independently defined and each has one 
and only one solution in the corresponding Sobolev space. The problem (13.121) is a 
standard second order elliptic equation for which many existing numerical methods 
can be applied. But the problem (13.211) requires a study in numerical methods, which 
is the central topic of this paper. 


4. The div-curl System with Tangential Boundary Condition. The goal 
of this section is to derive a variational formulation for the div-curl system with the 
tangential boundary condition dm. Recall that the complete problem reads: given 
g G with V ■ g = 0 in n, find the vector-valued function u such that 


(4.1) 


V • (mu) = /, 
V X u = g, 

u X n = X, 

(mu • Ui, l)r, = Pi, 


in n, 
in n, 
on r. 


The problem (14.11) can be interpreted as a constraint minimization problem with 
an object function given by 

= ^l|Vxv-g||g, 

subject to the following constraints 

r V • (mv) = /, in n, 

(4.2) <1 V X n = X, on r, 

[ (mv ■ YVi, l)ri = Pi, i = l,...,m. 

By introducing a Lagrange multiplier p, the corresponding variational problem seeks 
u G iL(curl; Lt) DH{div fl) and p G Lp{Ll) with u x n = x ou T and (mu-u^, l)ri = Pi 
for i = 1,..., m such that 

(Vxu, VxC) + (V • (mC),p) = (Vxg,C), 

(V • (mu),u>) = {f,w), 
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(4.3) 


VC e Y^(r!)niL(div^;L!), 
y wGL^{VI). 


















The problem (14.3p is the desired variational form for the div-curl system with 
tangential boundary condition dm. The structure of this variational problem is 
essentially the same as that of (13.211) arising from the normal boundary condition. 

Lemma 4.1. The solution to the variational problem 'Is unique for any 

connected domain fl. 


Proof. It suffices to show that all the solutions corresponding to the one with 
homogeneous data are trivial. To this end, let (u;p) be a solution of (14.3p with 
homogeneous data, then 


(4.4) 

(4.5) 

(4.6) 

(4.7) 


Vxu = 0, 
V • (^u) = 0, 

(^u • Ui, l)r, = 0, 
u X n = 0, 


in n, 
in n, 

z = 1,..., TO, 
on dn. 


From (|4.5I) and (|4.6L the conditions of Theorem 3.4 in m are satisfied for the vector¬ 
valued function /iu. Thus, there exists a vector potential 4> € [i7^(0)]^ such that 
fiu — Vxcf). It follows that 


(/ru,u) = (Vxc/),u) 

= (</), Vxu) -I- (n X </), u)r 
= n)r 

= 0 , 


by integration by parts 

by (14.41) and triple product property 

by dm 


which implies u = 0. Consequently, one has the following equation 


(V • (K),p) = 0, VC G Y^(fl) n i7(div^; fl), 

which leads to p = 0 by selecting a particular vector field C G Y^(fl) (~l 7L(div^;n) 
such that V- (pC) = P- Such a vector field is given by C = Vw, where w is the solution 
of the following equation 


V-(/rVw)=p, w|ro = 0, w\Ti = i = I, ■ ■ ■ ,'m. 

Here, { 7 ^} is a set of real numbers which can be tuned so that (/iVw ■ n^, l)ri = 0. □ 


5. A Model Problem. It is readily seen from Subsections 13.21 and Section 0] 
that the core problem of study for the div-curl system with either the normal or 
the tangential boundary conditions is one structured in the form of (13.211) and (14.31) . 
Thus, we shall consider a general problem in the form of (13.211) and (14.31) as follows. 
Assume the following data are given: g G [L^(r2)]^, / G T^(fl), C G [iV“5(r)]^, a 
symmetric and positive definite matrix k = (Kij(x))^^^ in the domain fl, and a set of 
real numbers fdi,i = 1 ,..., to. 


Problem 1. Find u G i7(curl; fl) (h i7(div^; fl) and p G L^{fl) such that 


(5.1) 


(kVxm, Vxu) + (V • (/rt)),p) = (g, v), 

(V • (Aiii),w;) = (/,u;), 

{p.u-n,,l)r, = |3^, 
ux n = Cl 


V w G Y^(H) (h i7(div^; fl), 

V w; G L'^ifl), 

Z = 1, . . . , TO, 
on r. 
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Here n is the unit outward normal direction to the boundary F. 

The constraint of {fiv ■ ni,l)ri = 0 for the test space Y^(r2) fl iF(div^;r2) is 
cumbersome in the design of numerical methods for the problem (EU- One possible 
remedy is to relax this constraint by using a Lagrange multiplier denoted by A = 
(Ai,..., Am)- The corresponding weak formulation seeks u e iF(curl; 0)n7L(div^; O), 
p G and A S R™ such that u x n = ^ on T and 

{ m 

(/cVxu,Vxv) + (V- {p^r),p) + = (g,v), 

(V • {pu),w) +'^{pu- n,, Si)r, = (/, w) + ^ /3iS,, 

for all V G iFo(curl;0) n 7F(div^;0), w G T^(0), and s G M’”. It is not hard to see 
that the problems (Ol) and (ED are equivalent to each other. 


For simplicity of analysis, throughout the paper, we assume that k and p are 
piecewise constant, symmetric and positive definite matrices on the domain O with 
respect to any finite element partitions to be specified in forthcoming sections. 

6. Weak Differential Operators. The model problem (15.11) is formulated with 
two principle differential operators: divergence and curl. This section shall introduce 
the notion of weak divergence operator for vector-valued functions of the form /rv. 
For completeness, we also review the definition for the weak curl operator. These 
weak differential operators shall be discretized by using polynomials, which leads to 
discretizations for the model problem ED- 

Let K G it he any open bounded domain with boundary dK . Denote by n the 
unit outward normal direction on dK. Let the space of weak vector-valued functions 
in K be given by 

ViK) = {v = {vo,v4 : VO G [L^iK)f, v„ G [L^idK)]^}, 

where Vg represents the value of v in the interior of K, and Vf, represents certain 
information of v on the boundary dK . There are two pieces of information of v on 
dK that are necessary for defining the variational formulation (15.11) : (1) the tangential 
component n x (v x n), and (2) the normal component of /iv on dK. The normal 
component of the vector /rv is given by (/rv • n)n. Intuitively speaking, the vector Vf, 
should carry those two pieces of orthogonal information by summing up (/iv • n)n and 
n X (v X n): 


Vb = (/rv • n)n -|- n x (v x n). 

It is easy to check the following identity: 

(6.1) Vb X n = (n X (v X n)) x n = v x n. 

6.1. Weak divergence. Following [21], for any v G V{K), we define the weak 
divergence of /iv, denoted by • (/iv), as a bounded linear functional in the Sobolev 
space H^{K) such that 

(Vu, ■ {pv),(p)K = -(/rvo, Vip)k + (vb ■ n, (p)dK, 
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V V? G H\K). 


The discrete weak divergence of /iv, denoted by ■ (/^v), is defined as the unique 

polynomial in satisfying 

(6.2) • (Mv),¥>)ic =-(Atvo, + (Vf, • n, (p e Pr{K). 

Assume that Vg is sufficiently smooth such that V • (fJ-vo) G L'^{K). By applying 
the integration by parts to the first term on the right-hand side of (I6.2I1 . we have 

(6.3) {yw,T,K ■ {f^'v),(p)K = (V • {fj.'Vo),(p)K + ((Vb - fiVo) ■ n, ip)aK, 
for any p G [PriK)]^- 

6.2. Weak curl. The weak curl of v G V{K) (see [2^), denoted by Vw x v, is 
defined as a bounded linear functional in the Sobolev space [H^{K)]^, such that 

(Vu, X V, p)k = (vg, V X p)k - (vb X n, p)dK, y pG 

The discrete weak curl of v G V{K), denoted by Vw,r,K x v, is defined as the 
unique polynomial in [Pr{K)]^, satisfying 

(6.4) (Vu,,r./f X V, p)k = (vg, V X p)k - (Vb X n, p)dK, V(^ G [Pr{K)]^. 

For sufficiently smooth vg such that V x vg G [L'^{K)]^, by applying the integra¬ 
tion by parts to the first term on the right-hand side of (16.4L we obtain 

(6.5) {yw,r,K X V, p)k = (V X Vg, p)k “ ((Vb - Vg) X tl, p)dK, 
for any p G [Pr{K)Y’ . 

7. Weak Galerkin Discretizations. A polyhedral partition of is a family 
of polyhedra {Tj : j = 1, 2,... } such that the two following conditions are satisfied: 
(1) the union of all polyhedra Tj is the domain and (2) the intersection of any 
two polyhedra, Tj n (i ^ j), is either empty or a common face of Tj and T^. 
Each partition cell Tj is called an element. A polyhedral partition with finite number 
of elements is called a finite element partition of n. Assume that {Tj'\j=i^,,,^N is 
a finite element partition of n that is shape-regular according to [30]. Denote by 
hr = diam(T) the diameter of the cell/element T, and h — maxT hr the meshsize 
of the partition Th = {Tj}j=i^...^N- Denote by £h tbe set of all faces in Th so that 
each e G fb is either on the boundary of D or shared by two elements. Denote by 
= £h\ d£l the set of all interior faces in £h. By definition, for each interior face 
e G £^, there are two elements Tj and i j, such that e = Tj r\Ti. 

Let A: > 1 be any integer. For each element T G Th, define the local finite element 
space as 

V{k,T) = {v = {vg,Vb} : Vg G [Pfe(T)]3,Vb G [Pfe(e)]3, e G {dTn£h)}. 

The global weak finite element space for the vector-component is given by 
(7.1) Vb = {v = {vg,vb} : v|t G V(fc,T), VblaTine = VblaTsne, T G 7b,e G f°}. 
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where 'Vb\dTjne is the value of on the face e as seen from the element Tj, j = 1,2. 
The finite element space for the Lagrange multiplier p is defined as 

Wh = {q-. q€ L^in), q\T S Pk-i{T),T € %}■ 

The discrete weak divergence and the discrete weak curl (Vu,_fc_ix) 

can be computed by using (j6.2l) and (16.41) on each element; i.e., 

• (/rv))|7’ =Vu,,fe-i,T • (mv|t), V G V/j, 

(Vu,,fe-i X v)|t =Vu,,fe-i,T X (v|t), V G V?,. 

For simplicity of notation, we shall drop the subscript k — 1 from the notations 
and (Vu;,fe-ix) from now on. 

Introduce the following bilinear forms 

(7.2) a(v, w) =(/cVu, X V, Vu, X w)h + s(v, w), 

(7.3) b{-v,q) ■ {p'v),q)h, 

where 

(kVu, X V, Vu, X vf)h = ^ (kV^ X V, Vu, X w)t, 

T&Th 

(Vu, • ipv),q)h = ^ (V„ • {pv),q)T, 

Ten 

(7.4) s(v, w) = ^ /iy^((^vo - vf,) • n, (^wq - w;,) • n)aT 

Ten 

+ ^ /it^((vo - Vfe) X n, (wo - Wf,) x n)aT. 

Ten 


We are now in a position to describe a finite element method for the model problem 
given in dSU) and (15.211 . Consider first the variational formulation (15.211 . Observe that 
the test space is i7o(curl; fl) 0 i7(div^; fl). The corresponding analogue in the weak 
Galerkin setting is the following weak finite element space 

(7.5) = {v = {vo,Vb} G : v^xn^OonT}. 

Algorithm 1. (weak Galerkin for the problem (15.21) 1 Find Uh = {?xo: G V^, 

Ph S Wh, and X G with Ub x n ~ Qb^ on T such that 

m 

(7.6) a{uh,v)+h{v,ph) + ^{vb-ni,Xi)Ti = {^,VQ), V ?; = {vq, Vb} G V((, 

2 = 1 

m m 

(7.7) b{uh,w)+^{ub-ni,Si)Ti = {f,w)+'^l3iSi, V n; G IFb, s G M™. 

i —1 2=1 

Here Qb^ is the usual L'^-projection to [Pfc(e)]^ for each boundary face e G (T 0 £h)- 

The Lagrange multiplier A can be eliminated from the formulation dZUl-dLll) if 
the test space V° is replaced by 

(7.8) = {v = {vo, Vb} G : (vb • n^, l)r, =0, i = 1,..., m}. 
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The following is the corresponding weak Galerkin finite element scheme. 

Algorithm 2. (weak Galerkin for the problem (j5.ll) ') Find Uh = {tio,ui,} G 
and ph G Wh with Uf, x n — on T such that 


(7.9) 

(7.10) 

(7.11) 


a{uh,v) + h{v,ph) = (g, i>o), 
b{uh,w) = {f,w), 

{ub ■ Ui, l)ri = /3i, 


V 11 = {wciib} e Uj), 
y w G Wh, 
i = 1,... ,m. 


Note that the Algorithms [T] and [2] are equivalent in the sense that both give the 
same numerical solution and ph- Thus, it is sufficient to develop a convergence 
theory for Algorithm [2] only. 


8. Existence and Uniqueness. The goal of this section is to show that the 
WG Algorithm has one and only one solution. 

We first introduce a topology in the weak finite element space by defining a 
semi-norm as follows 


( 8 . 1 ) 


T€Th 

1 

+ h^^ll/ivo • n - Vb • n||ay -I- /i;^^||vo x n - Vb x njlg^^ ^ 


For convenience, we set 

|||v|f := a(v,v) = ^ (kV^ 


( 8 . 2 ) 


TgTh 


hj,'^\\pvQ 


X V, Vu, X v)t 

• n - Vb • n||ay -I- ||vo x n 



1 . 

2 


In the finite element space Wh, we introduce a mesh-dependent norm 

m 

(8-3) WlWwh — '^2 ll^gllr + ^ lllgielle + l|g ~ QiWri^ 

TeTh e<^e° i=0 

where [gfg stands for the jump of q on the interior face e G £^, go = 0, and qi is the 
average of q on the connected boundary component Fj, i = 1,..., m. 

Lemma 8.1. Assume that the domain fl is eonnected. Then, the semi-norm ||| • |||j^ 
defined as in H8.1\} defines a norm in the linear space U®. 

Proof. It suffices to verify the positivity property for ||| ■ ||| ^. To this end, assume 
that |||v|||^ = 0 for some v G U°. Thus, 

(8.4) Vm X V = 0, in T, 

(8.5) Vuj • {pv) = 0, in T, 

(8.6) pvo ■ n — Vb ■ n = 0, on dT 

(8.7) Vq X n — Vb X n = 0, on dT. 
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Using (EH), (EH) and EH, for any ip e [Pfe_i(T)]^, we have 


0 =(Vu, X V, (/j)t 

=(V X Vo, p)t - {{yb - Vo) X n, p)aT 
=(V X vo,(/3)t- 


It follows that V X Vo = 0 on each element T G Th, which along with (18.71) implies 
Vo G i7(curl;f2) and 


( 8 . 8 ) 


V X Vo = 0, in Q. 


Next, using (18.51) . (18.61) and (16.31) . for any p G we have 


0 =( Vu , • {p'v),p)t 

= (V • {p'Vo),p)t + ((vf, - /rvo) • n, p)9T 
= (V • {pvo),p)t. 


It follows that V ■ (^vo) = 0 on each element T gTh, which, together with (18.6|) . gives 
rise to ^vo G i7(div; fl) and 


(8.9) 


V • ifJ.'Vo) = 0, in fl. 


Combining (I8.6I) - (I8.7I) with the fact that v G yields 

Vo X n = 0, on F, 

(/rvo ■ tii, I)ri = 0, i = l,...,TO, 

which, along with (|8.8p and (18.9L implies Vo = 0 and hence v^ = 0; see the proof of 
Lemma Q] for details. □ 

Theorem 8.2. The weak Galerkin finite element algorithms\^ and\^have one 
and only one solution. 

Proof. Since the Algorithms [T] and [5] are equivalent, then it suffices to deal with 
Algorithms [2J As the number of unknowns is the same as the number of equations in 
(I7.9[) - (I7.11[) . the existence of solution is equivalent to the uniqueness. 


Let G Vh x Wh, j = 1,2, be two solutions of (I7.9I) - (I7.1I|) . and set 



It is clear that {zh;jh) G Uj) x Wh satisfies 


( 8 . 10 ) 

( 8 . 11 ) 


a{zh, v) + 6(v, 7 /t) = 0, V V = {vo, v4 G U°, 

b{zh, w) = 0, \/ w G Wh. 


By letting v = Zh in (|8.10l) and w = y/i in (|8.11l) . we obtain 

a{zh,zh)=0, Vu, • (/rz,i) = 0, 

which leads to |||zft,|||i = 0. It follows from Lemma f8.ll that Zh = 0, and thus = 
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To show 7 /j = 0, we use (18.101) and z^j = 0 to obtain 


b(v,'yh) = 0, V V = {vo.vb} G U^. 

From the definition of weak divergence (IQ) . 

b(^,lh)= ^ (V™ • (^v),7,i)t 
T&Th 

(8.12) = -(^vo, V7 ;j)t + (v6 • n,7?,)aT 

TeTh 

T^Th e.^Sh 

where is a prescribed orientation of e. By letting v = {—h^V^h', h5e\"ih\T^e} 
with i5e = 0 when e C Fi, i = and (5e = 1 otherwise, we see that v G U°. 

Substituting this into (18.121) yields 

Q = h{'v,-fh) = h^ ^{^JS/-|h,S/-fh)T+ h ^ Ill 7 /illl^ 

T&Th eefguro 

It follows that 7 /i = 0, and thus This completes the proof of uniqueness. 

□ 


9. Error Equations. Let Qo be the projection onto [Pk{T)]^, T GTh, and 
Qb the projection onto [Pk{e)]^, e G dT l~l Eh- Denote by Qh the projection 
onto the weak finite element space V/i such that on each element T gTh, 

(9.1) (Q/iU)|t = {QqU, Qfeu}, 

where 

(9.2) QfeU = Qb(iiu • n)n + (2{,(n X (u X n)). 

Note that n x (u x n) = u — (u • n)n is the tangential component of the vector u 
on the boundary of the element. When p = I is the identity matrix, {pu ■ n)n is the 
normal component of u. In general, (^tu • n)n + n x (u x n) is not a decomposition 
of the vector u restricted on dT. 

Denote by Qh and the projections onto Pk-i{T) and [Pk-i{T)]^, respec¬ 
tively. 

Lemma 9.1. \2(A \31f The projection operators Qh, Qh, bind Qh satisfy the fol¬ 
lowing commutative identities: 

(9.3) yw-{pQhv)=Qhy-{pv), G iL(div^;fI), 

(9.4) X-{Qhv) ^Qhi^ v), t) G i?(curl; D). 


Proof It suffices to verify (lOD on each element T € Th- To this end, using the 
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definition (lOl) for the discrete weak divergence and (ESI), we obtain 

{Vyj ■ {fiQhv),(p)T = -(/rQoV, V(p)t + ((Qbv) • n, (^)ar 
= -(mQoV, Vif)T + {Qb{^J■^v ■ n), ip)oT 
= -{nv, Vip)t + ifJ-'v ■ n, (p)dT 
= (V • {fiv),if )T 
= (Q/iV • {fj.v),(p)T 

for all (fi G Pk-i{T). Thus, the identity (19.3|) holds true. A similar argument can be 
applied to verify (19.41) . □ 

Let {uh',Ph) = ({uo, ui,}\ph) G 'Vfi X Wh be the WG finite element solution arising 
from Algorithm [21 and (u;p) be the solution of the continuous problem (15.11) or (15.21) . 
The error functions are given by 

(9.5) Bh = {eo, Bf,} = {Qou - Uq, - U;,}, 

(9.6) eh = QhP-Ph- 

Lemma 9.2. Assume that (w\p) G iL(curl;n) x L^(n) is sufficiently smooth on 
each element T gTh satisfying 

(9.7) V X (kV X to) — pS7p = p, in fl, 

(9.8) pIto = 0, pIt; = const, i = 1,... ,m. 

Denote by Qhp the projection of p in the finite element space Wh- Then, 

(9.9) (/cVu, X (Qhw), Vyj X v)h + (Vu, • (pv), Qhp)h = {p, I’d) + lw{v) + 9p{v), 

for all V G U((. Here l-w{v) and 0p{v) are two functionals in the linear space given 
by 

(9.10) l.u,{v) = '^{{Qh- -f)(KV X to), (vo - Vh) X n}9T, 

Ten 

(9.11) 0p(v) = ^ (p - Qhp, (pvQ - Vb) ■ n} 9 T- 

Ten 

Proof. It follows from (16.51) with (p = kV^j x {Qh'W') that 
(Vu,xv, kVu, X {QhVf))T = 

(V X Vo, kVu, X - ((v& - vq) x n, kV^, x {Qh'^))dT- 

Using dSAl), the above equation can be rewritten as 
(KVu,x((5?,v^f), Vu, X v)t = 

(kV X w, V X vo)r + (Q/t(KV x w), (vq - V(,) x n)oT- 
Applying the integration by parts to the first term on the right-hand side yields 
(/cVu,x((5hw), X v)t 

=(Vx(kVxw), Vo)t — (kVxw,Vo x n)gr 

+ (Q/i(«V X w), (vo - Vb) X n)aT 
= (Vx(kVxw), vo)t — {KVxw,Vb X n)aT 
+ ((Q/i - -f)(KV X w), (vq - Vb) X n)aT- 
16 


(9.12) 


Using (lOD with ip = QhP and the usual integration by parts, we obtain 
(V^ • (/rv), Qhp)T 

= (V • (/rvo), Qhp)T + ((vf, - ^vo) • n, Qhp)dT 
(9.13) = (V • (pvo), p)t + {(vb - p-vo) • n, Qhp)dT 

= - {pvo, Vp)t + (/rvo • n, p)oT + ((vf, - /rvo) ■ n, Qhp)dT 
= - (vo, pV p)t + ((vb - p^fo) ■ n, QhP - p)dT + (vb • n, p)dT- 


Summing (19.121) over all the elements T GTh yields 

{K.VyjX{QhVf),Vyj X v)b = ^ (Vx(kVxw), vo)t 


(9.14) 


T(^Th 

+ E m, - I){kV X w), (vo - Vb) X n)QT, 
T&n 


where we have used two properties: (1) the cancelation property for the boundary 
integrals on interior faces, and (2) the fact that Vb x n = 0 on F. Similarly, summing 
(19.131) over all the elements T gTh, obtain 


(9.15) 


(V^ • (pv), Qhp)h = - (vo, pyp) + ((mvo - Vb) ■ n, p - Qhp)dT 

TGTh 

+ (v&'n,p)e. 

eeShnr 


The third term on the right-hand side of (19.151) vanishes if v S U° and p satisfies the 
boundary condition (19.81) . Thus, the equation (19.91) holds true by adding up (19.141) 
and (|9.15l) . This completes the proof of the lemma. 0 

Theorem 9.3. Let {u]p) be the solution of the model problem \5.1\) or i5.2\} and 
(uh]Ph) be its numerical solution arising from the WG finite element scheme 
([7.11| ). Let the error functions eh and Ch be defined by l,9.5l) - l,9.dl) . Then, eh € 11° 
and the following error equations hold true 

(9.16) a{eh,v)+b{v,eh) = (Pu,p{v), Vi; e U((, 

(9.17) b{eh,q)=0, yqeWh, 

where 


(9.18) 


Pu,p{v) = lu{v) + 9p{v) + s{QhU, v). 


Proof. Consider the model problem (15.21) . as m IS equivalent to (|h,^j) . Let 
(u;p; A) be the solution of this model problem. It is not hard to see that the following 
holds true: 

V X (kV X u) — pVp = g, in Q, 
pIeo = 0, p|ri =-Ai, i = l,...,m. 

Thus, by Lemma 19.21 we have 

(kVu, X (Qbu), Vyj X v)b -f (Vu, • (pv), Qhp)h = (g, Vo) -f Zu(v) -f 9p{v) 


17 














for all V e U°. It follows that 

(9.19) a{QhU, v) + 6(v, Qhp) = (g, vq) + /u(v) + 9p{v) + s{QhU, v). 

Subtracting (17.91) from (I9.19P gives the first error equation (19.161) . 

Next, from the second equation in (15.11) and the commutative relation (19.3L 

(9.20) if, (?) = (V • (/ru), q)h = {Qh^ ■ (mu), q)h = (V^ • {pQh^), q)h- 

The difference of (19.201) and (17.101) yields the second error equation (I9.17p . This 
completes the proof. 0 


10. The inf-sup Condition. For any q S Wh, define a finite element function 
V, e by V, = hv,,;,}: 


( 10 . 1 ) 


f [ 9 I Ue, on e e 

\ {q- qi) n*, on e G n T*, i = 0,1,..., m. 


Recall that [(?[ is the jump of q on the corresponding face e G ne is a prescribed 
orientation of e, n^ is the outward normal direction on the connected component T^, 
go = 0, and qi is the average of g on T^, i = 1,... ,m. 

For any v = {vq; Vf,} G U° , from the definition of weak divergence ()6.2I) . we have 


^(v,g) = 

TdTh 

= X! -(mvo, Vg)r + (v& • n,g)aT 

TdTh 

m 

= - X! (mvo, Vg)T Ue, [gl)e + ' Uz,(?)ri- 

TGTh ee£g i=0 


Note that (v;, • n^, l)ri = 0 for ? = 1,..., m. Thus, 

m 

(10.2) 6(v,g) = - (?xvo, Vg)T + Y + ' “^’9“ *)ri- 

T^Th ee£° i=o 


Lemma 10.1. (inf-sup condition) For any q G Wh, there exists a finite element 
function Vg G U° such that 

m 

(10.3) b{vg,q)=h^ XI + ^ + 

TGTh e&£° i=0 

(10-4) |||?;9||| < llglk^. 


Proof. For any g G Wh, define Vg^t by (110.11) and set Vg = {—/i^Vg; hvg^b}- On 
the boundary F, the vector Vg^t, is parallel to the normal director n. Thus, we have 
"Vq^b X n = 0 on F. Moreover, on each connected component F^, we have 


{vq^b ■ rii, l)ri = J {q-q^) = 0, i = 1,2,...,', 
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Thus, v,j e U°. 

Now, by taking v = in p0.2|l . we obtain 


m 

b{vq,q) = h^ (mVq, Vg)T + h ^ IIfcl lie + /»XI Il9“9i|lr,> 

TeTh e(^S° i=0 


which verifies the identity (110.31) . 

To derive (I10.4L we consider the following decomposition 


fq = 


.( 2 ) 


where = —{h^Vq] 0} and = {0; It suffices to establish (110.41) for 

and independently. 

From the semi-norm definition (18.2L we have 

lllv^^^lf = X ^ X vW)g, 

(10.5) TdTh 


+ hj-^Wli^^Vq ■ ri\\^QT + /iy^||h^Vg x nUgy. 

The definition (16.41) for the discrete weak curl implies 

(V^ xv«,(^)t = -/i2(V<7,Vxvj)t, V(^e [Pfe_i(r)]3. 

It follows from the inverse inequality that 

l|V» xvW||T<h||Vg||T. 

Substituting the above into (110.51) and then using the trace inequality (111.31) yields 

</i"||Vg|||, 

which verifies the estimate (110.41) for Vg^^ 

For Vg^^, we again use the semi-norm definition (18.21) to obtain 

(10.6) |||v(^)|f = X («;Vu,xv(^\VujXv(^))T + /iy^||/iVg,{,-n||aT-fhy^||/iVg,(,xn||ay. 

TdTh 

Since Vg ^ is parallel to n, then Vg_{, x n = 0 on dT. In addition, the definition (16.41) 

f 2) 

for the discrete weak curl implies Vu, x Vg = 0 as 

(Vu, X w^q\q})T = (0,Vxip)t - h{vq^b X n,(f)9T = 0, V £ [Pk-i{T)f. 

Thus, it follows from (110.61) and (110.11) that 


ftllr-i 


lllvff X Ill9llle + Xll9-' 

\ee£« i=0 

which verifies the estimate (110.41) for Vg . This completes the proof of the lemma. □ 
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11. Error Analysis. Based on the error equations shown as in Theorem 19.31 
and the inf-sup condition in the previous section , we shall derive an estimate for the 
error terms eh and eh taken as the difference of the WG finite element solution and 
the projection of the exact solution. 

11.1. Some technical inequalities. Assume that the finite element partition 
of fl is shape regular as defined in [30]. Let T S 7h be an element with e as a face. 
The trace inequality holds true: 

( 11 . 1 ) mi<ihT^urT+hT\mrT), 

If (ji is a polynomial, the inverse inequality holds true: 

(11-2) l|V(/>||T</l?i0||T. 

From (111.11) and (I11.2L we have 

(11-3) 

Lemma 11.1. JWf Let k > 1 be the order of the WG finite elements, and 1 < 
r < k. Let w € p G and 0 < m < 1. There holds 


(11.4) 

(11.5) 

( 11 . 6 ) 




2 

r+1) 


Ten 


4™I|V XW- Qh(y X w)\\l^ < 

Ten 

E hl^\\p-QhPrn^<h^''\\p\\l. 

Ten 


In the WG finite element space V/j, we introduce a semi-norm as follows 

(11.7) |v|i,,,= I ^ /i^^ll(vo-vb) xn||^r + ^?^ll(MVo-V6)-n|||7.| . 

XTen / 

Lemma 11.2. Assume that the finite element partition Tf of LI is shape regular 
as defined in \3Uf and 1 < r < k. Let w S [iL'’“'"^(r2)]^ and p G Then, we have 

(11-8) \s{QhW,v)\ < h^\\w\\r+l \v\i^h, 

(11-9) \lw{v)\ < h''\\w\\r+l \v\i^h, 

(11-10) \Sp{v)\ < h''\\p\\r \v\i^h, 

for any v G V/j. Here lw{-) and 9p{-) are defined in i9.10\) and 

Proof. Recall from (17.41) that the stability term can be decomposed into two parts: 



s(v, 

, w) = Si(v 

■,w) -h S 2 (v,w), 

where 




(11.11) 

Si(v,w) = 

(('^O 

Vfe) X n, (wo - Wfe) X n)oT, 


Ten 



(11.12) 

S 2 (v,w) = Y 

hT^iip^o - 

- Vfe) • n, {pwo - Wh) ■ n)9T 


Ten 
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To prove (lll.8|) . it suffices to derive the estimate (111.81) for Si(-, •) and S 2 (-, •) sepa¬ 
rately. To this end, we use the Cauchy-Schwarz inequality, the trace inequality (lll.ip 
and the estimate (111.41) to obtain 

|si(Q/tW,v)| = I ^ h~^{QoVf X n- (^^(w x n), (vq - v;,) x 
T eTh 

= I X! ^t^(Qow X n - w X n, (vo - vf,) X n)aT 
TGTh 

^t^Qow-/ iT^II(vo - Vf,) xn|||-r)" 

TGTh TGTh 

- ( X! ^t^IIQow - will -h IQow - wIIj.) " |v|i,?, 

T&Th 

< /l’'||w||r+l |v|i,?,. 

To derive ()11.8p for S2(-,-)) again use the Cauchy-Schwarz inequality, the trace 
inequality (111.11) and the estimate (111.41) to obtain 

|s 2 (Q/iW,v)| = I Y, • n - (3h(/rw • n), (/rvo - Vh) • n)aT 

TGTfc 

= I X! • n - /iw • n, (^vq - Vf,) • n)aT 

TGTfc 

^ ( X! w)|||t) ' ( Y ^Ti(AiVo-Vf,)-n|||r)" 

TGTh TgTh 

- ( X! ^t^II<3oW- w|||-H IQow-w||-r) " |v|i,,i 

TGTh 

< /l’'||w||r+l |v|l,/t. 


As to (I11.9L we use the Cauchy-Schwarz inequality, the trace inequality l|ll.ip 
and the estimate (111.51) to obtain 


l^w(v)| = 


Y - d)(KV X w), (vo - Vfe) X n))oT 


TgTh 


< 


( Y ^T\\{Qh-I)(y xw)||ot)"( Y ^t^II(vo-V 6) xnlllr) 


TgTh 

</l’'||w||r+l \v\i^h- 


TGTh 


Finally, we use the Cauchy-Schwarz inequality, the trace inequality (jll.ip and the 
estimate (lll.bp to obtain 


|6'p(v)| = 


Y ^p~ (^'"0 “ '"f’) ■ 


dT 


TGTh 


< 


( Y ^r||p- Q^^pIIot) " ( Y ^t^IK/^vo - Vfa) • nlllr) 

TgTk TgTh 
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This completes the proof of the lemma. □ 

11.2. Error estimates. Recall that ||| • |||j^ defines a norm in the finite element 
space U^. This norm can be regarded as a discrete i?o(curl) fl i7(div^)-norm under 
which the error function shall be measured. 

Theorem 11.3. Assume that k > 1 be the order of the WG finite element scheme 
vm-urm- Let iu;p) G x H^{n) be the solution of the model problem 

fsnp and {uh'jPh) be the WG finite element solution arising from l[7.9^ - /i7.11\ ). Then, 
we have 

(11.13) WQhU- Uhii + WQhP - PhWwh ^ ^^(ll'“IU+i + IblU). 


Proof. From Theorem 19.31 the error functions et = Qh^ — u/i and et = QhP — Ph 
satisfy the equations (I9.16I) - (I9.17I) . By letting v = in l|9.16|) and then using (j9.17|) . 
we obtain 

(11-14) a{eh,eh) = (pu,p{eh)- 

The right-hand side can be estimated by using Lemma 111.21 as follows 

|<7’u,p(eh)| < /i''(||u||fe+i -b |b|U)|e,i|i,,,. 

Substituting the above into (111.141) yields 

a{eh,eh) < /i'‘(||u||fc+i -b |b||fc)|e;i|i,?,, 

which, together with \eh\f < a(e/j,e/i), leads to 

(11-15) aien,e^,)^/^ <h\\\u\\k+l + \\p\\k)■ 

lSlext, it follows from the equation (19.171) that • {peh) = 0. Thus, 

ll|e/i||li = |||eft||| < a{eh,ehy^^ 

Combining the above inequality with (111.151) gives 

(11-16) ll|e/i||li < h^i\\u\\k+i + Ibllfe)- 

The error function Ch can be estimated by using the inf-sup condition derived in 
Lemma riO.il To this end, from the equation (19.161) . we have 

(11-17) b{v,€h) = Pu,p{v) - a{eh,v). 

By letting v = in (111.171) , we arrive at 

l|e/«llk - l7^u.p(v,J| + |a(e,„v,J|. 

Using Lemma [11.21 and the error estimate (Ill.lSp we obtain 

lk.llk^^''(ll’^IU+i + lblU)IKJII, 
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which, together with (I10.4L leads to 

W^hWw^ < h'=(||u||fc+i + ||p||fe). 

This completes the proof of the theorem. □ 

The mesh-dependent norm UgllvK;, in the finite element space Wh is a scaled dis¬ 
crete norm for piecewise smooth functions. For the lowest order WG element (i.e., 
piecewise linear for u and piecewise constant for p), the error estimate in Theorem 
lll.3l does not give any convergence for the approximation of p. However, it is possible 
to replace the || • \\wh norm by the standard norm in the error estimate (111.131) if 
the solution of the following second order elliptic problem is i7^-regular: 

-V • (/xV$) = e/i, inH, 

= Cti, i = 0, 1, ... , 771 , 

where ao = 0 and is a set of constants such that (V$ • rii, l)ri = 0 for i = 1,..., tti. 
In fact, for the solution of (I11.18L it can be seen that V$ S n H(div^;H). 

Moreover, the projection z = (see (19.11) for its definition) is a finite element 

function in U)(. The desired error estimate for in LF’iVl) can then be obtained by 
taking v = z in the error equation (19.161) . Details are left to interested readers as an 
exercise. 


11.3. -error estimates. To derive an error estimate for e^,, we consider 
the dual problem of seeking i/) € iJo(curl; H) fl i?(div^; H) and r G such that 

(kVxt/j, Vxv)-h (V • (/iv),T) = (eo,v), V v G Y^(H) n i7(div^; H), 

(11.19) (V • (/7i/)), w) = 0, V Tc G T2(D), 

{pip ■ rii, l)ri = 0 , 7 = 1 ,... ,m. 

Assume that the dual problem (111.191) has the [iJ^(r2)]^ x i?^(f2)-regularity property 
in the sense that the solution G [i?^(D)]^ x and satisfies the following a 

priori estimate: 

(11-20) llV’lb + ||r||i < lleollo. 


By using a Lagrange multiplier 7 = ( 71 ,..., 7 ^) G M’”, the dual problem (I11.19|) 
can be rewritten in an equivalent form as follows. Find xp G iLo(curl; H) niL(div^; H), 
T G L^(D), and 7 G K.™ such that 


( 11 . 21 ) 


(kVx-0, Vxv) -h (V ■ {pv),T) + ^{pv ■ ni,7i)r, = (eo,v), 

2=1 

m 

(V • {pxp),w) + '^{pxp- rii, Si)r, = 0, 


for all V G iLo(curl; H) n H{divp H), ic G T^(H), and s G R™. 

Theorem 11.4. Let k > 1 be the order of the WG scheme [m-r7G7\)- Let 
{u;p) G [iL^+^(H)]^ X iL^(H) and {uh',Ph) G 'Vh x Wh be the solutions of the problem 
m\> and respectively. Then, the following estimate holds true 


( 11 . 22 ) 
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Proof. From the first equation of (lll.21|) . we see that the the equation (I9.7|l is 
satisfied by (w; p) = r) with ry = Gq. In addition, the boundary condition (19. 8p is 
verified by 


'Tiro = 0, r|r, = -7i, f = 1,..., m. 

Thus, by using (19.91) in Lemma [Q] with v — G we obtain 


(11.23) 


Iqu - uqIP =(/cVu, X {Qhip), Vu, X eh)h 

F (Vuj • (y/e/j), QhT'jh, ^Tpiph)- 


Note that the error equation (I9.17P implies V^, ■ [peh) = 0 so that (I11.23P can be 
simplified as 

(11-24) IIQoU - uolP = a{Qhip, e^) - (^v.T(e?i), 

where 

Pip — ^riph) F ^ipi^^h) F s((5/i'0, 

From (19.31) and the second equation of (111.191) . we have 

w ' {pQh'^f^h'jh 

= {Qh^ ■ {p4’),<^h)h = 0, 

which, together with (19.161) and (I11.24F leads to 

IIQoU - uqIP = a{eh,Qh'if) F h{Qh'if, eh) - PipA^h) 


(11.25) 


— Pu,p{Q h'^) Pip ,r(.^h) • 


The two terms on the righ-hand side of (I11.25P can be estimated as follows. First, 
by letting r = 1, v = e?i, and (w; p) = {xp; r) in Lemma lll.21 we obtain 

(11-26) \Pip,T{eh)\ < h{\\'ip \\2 F ||r||i) < h\\eo\\ |||e/,|||, 

where we have used (jll.20|) in the second inequality. Next, by letting r = k, v = QhXp, 
and (w; p) = (u;p) in Lemma lll.21 we arrive at 

(11-27) \Pii,p{Qh'4’)\ ^ ^^(l|u||fe+i F \\p\\k) \Qhip\i,h- 

To estimate \Qhip\i,h, we recall from pi.71) and the definition (19.11) of Qh that 
(11.28) IQh'ipfiM = ^T^\\iQo'>P - Qb^) X n||9T F hAWipQo'ip - Qh^) • 

TGTh 

By using (Ha), we have 


(Qh'*/’) X n = QbA X (•0 X n)) X n 
= QbXp X n, 

(Qb0) • n = Qb{pxp - n). 
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Thus, 


(Qo'0 - ^bi>) X n = ((5o'0 - Qbip) X n, 
- ^bip) ■ n = n{Qoi> - Qbip) ■ n. 


Substituting the above identities into (jll.281) yields 

IQb-iPll^h = K^WiQoi’ - Qb-ip) X n\\lj. + h-^WiiiQoxjj - Qb-ip) ■ n||^y 

TGTh 

< hrp^WQo'tp — QbipWoT 

T&Th 

< ^T^WQoi’-'fPWdT 

TeTh 

< Y, i^T^Wi’ - Qoi’Wr + l|V(^ - QoV’)IIt) 

TeTh 

< Y 4iivviii 

TeTh 

<h^^P\\l 

where we have used the property of Qb in line 3, the trace inequality (111.11) in line 
4, and the usual approximation property for the projection operator Qq in line 5. 
Inserting the above estimate into (jll.271) and then using (lll.20|) . we obtain 

(11-29) \^M.piQh'ip)\ < h'=+i(||u|| fe+i + IblU) ll®o||- 

Now, combining (111.251) with (I11.26P and (jll.291) . we arrive at 

IIQoU - uof < (h\\\eh\\\ + /i'=+i(||u||fe+i + ||p||fc)) ||eo||. 

It follows that 

IIQou - uoll < h|||e^||| + /i'=+^(||u|| fc+i + Ibllfc), 

which, together with Theorem 111.31 and the fact that |||ef,||| < |||e?,|||j^, completes the 
proof of the theorem. □ 

The WG finite element solution Uh = {uo;Uf,} consists of two components on 
each element T G Th' (1) the element interior value Uq, and (2) the element boundary 
value Uf). The rest of this section is devoted to an error analysis for Uf,. To this end, 
we introduce the following topology in the finite element space V/, 

r - 

IWbh^^i^Y^'^ V= {voiVb} e V?,. 

^ TeTfc 

It is clear that the above defines an L^-like norm for the face variable v;,. 

Theorem 11.5. Let k > 1 he the order of the WG scheme Let 

{u]p) € X H^{Vl) and {uh]Ph) G V/j x Wh be the solution of the problem 

and CT-fCTP, respectively. Then, we have 

W^bU-UbWe,, < /i'=+^(||n||fe+i + ||p||fe). 


(11.30) 
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Here Qf,tt is the projection of u on each face e G £h given by 

Proof Note that — u;, = e;, is the error of the WG finite element solution on 
£h - the set of all element faces in 7^. It follows from the triangle inequality that on 
each dT 

|ef,p = |ef, • np + |e& x np 

< 2\{peo - eb) ■ np + 2|(eo - e^) x np 
+ 2\fieo • n|^ + 2|eo x n|^. 

Summing over all the elements and then using (ED yields 

[ |ebpds< |||e/j|||^ + ^ f (|^eo • np + |eo x np)ds 
TGTh TGTh 

<||K|||^ + h-2||eof, 

where we have used the trace inequality (jll.3|) in the second line. The last inequality 
further leads to the following estimate: 

/ |e&pds </i^|||e,i|||^ + ||eo|p, 

TeTh 

which, together with the error estimates (111.131) and (111.2211 . implies the desired 
estimate (ursQi). □ 
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